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Abstract 

We calculate an s-wave amplitude matrix for all the possible 2-to-2 body scalar boson elastic scatterings in 
models with three scalar doublets, including contributions from the longitudinal component of weak gauge 
bosons via the Equivalence Theorem Approximation. Specifically, we concentrate on the two cases with 
two[one] active plus one[two] inert doublet fields, referred to as I(1+2)HDM[I(2+1)HDM], under CP conser¬ 
vation. We obtain three analytically irreducible sub-matrices with the 3x3 form and eighteen eigenvalues 
for the amplitude matrix as an independent set, where the same formula can be applied to both models. By 
requiring a perturbative unitarity condition, we can constrain the magnitude of quartic coupling constants 
in the Higgs potential. This constraint, in particular in the I(1-|-2)HDM, can be translated into a bound on 
masses of extra active scalar bosons. Furthermore, when Standard Model-like Higgs boson couplings with 
weak gauge bosons are deviated from the Standard Model predictions, the unitarity condition provides an 
upper limit on the masses. We find that stronger upper bounds on the masses of the active CP-even and 
CP-odd Higgs bosons are obtained under the constraints from the unitarity and vacuum stability conditions, 
as well as the electroweak S, T and U parameters, as compared to those in 2-Higgs Doublet Models with a 
softly-broken Z 2 symmetry. 
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I. INTRODUCTION 


After the completion of the operations for LHC Run-I, the most important legacy for the 


Q. 


whole of particle physics is the existence of a Higgs boson with a mass of about 126 GeV 
Furthermore, as a result of combined data from various Higgs search channels, it was established 
that the observed Higgs boson has properties which are consistent with those of the Standard 
Model (SM) Higgs state Q]. 

However, it is natural to be left with a key question. Namely, what is the actual number of 
scalar doublet fields present in Nature? Although the Higgs sector in the SM is constructed of only 
one isospin doublet field, that is not supported by a fundamental principle. In fact, it is well known 
that multi-doublet models can be considered so that the Higgs boson phenomenology is consistent 
with current LHC data by arranging additional parameters suitably, and they are often introduced 
new physics models beyond the SM. 

For example. Supersymmetry-extended models require at least two doublets in their Higgs sector 


a. 


for gauge anomaly cancellation and to construct the Yukawa^Lagrangian Q] . A second doublet held 


fl, 


is introduced also in many radiative neutrino mass models . The Inert Doublet Model (IDM) 
aka a 2-Higgs Doublet Model (2HDM) with one active and one inert doublet, has been proposed 
to explain the existence of Dark Matter (DM), by depriving a Vacuum Expectation Value (VEV) 
of the second doublet field and making it non-interacting with SM fields except for the (thereby 
active) Higgs boson itself and the gauge states^. Eurthermore, an additional source of CP-violation 
is obtained in a Higgs sector with multrole doublets which is required to realize a successful scenario 
for Electro-Weak (EW) baryogenesis (6|. 

Properties of an enlarged Higgs sector clearly depend on the new physics scenario embedding 
it, so that the determination of the true Higgs boson dynamics is extremely important to extract 
information on the new physics model. In particular, if a multi-doublet structure is chosen for the 
Higgs sector, what is a mass scale of the additional active and/or inert (pseudo)scalar bosons? How 
can we extract information on the total number of doublet fields? Needless to say, the direct way to 
probe a multi-doublet structure is discovering extra (pseudo)scalar bosons at collider experiments. 
Por example, even in the simplest case of 2HDMs, there appear a pair of singly-charged a 

CP-odd (A) and a CP-even (H) Higgs boson in addition to the SM-like Higgs boson (h). Clearly, 


^ This can be simply realized by imposing an additional unbroken discrete Z2 symmetry into the Higgs sector, where 
all the SM particles [the second doublet] are assigned to be even [odd] under the Z2 symmetry. The particles with 
the odd parity are produced in pairs and the lightest one amongst these is stable, hence a DM candidate. 
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the most important property to find these extra scalar bosons is their masses. 

Although, in general, these masses are taken as free parameters, we may be able to constrain 
their corresponding parameter space by taking into account theoretical constraints: chiefly, pertur¬ 
bative unitarity and vacuum stability. In 2HDMs, the unitarity and vacuum stability bounds have 
been discussed in Refs. [?, 8| and in Ref. 9|, respectively. By using them, dimensionless param¬ 
eters in the scalar potential can be constrained, and it can be translated into limits on the mass 
parameters for the extra (pseudo)scalar bosons. In particular, when there is a non-zero mixing 
between h and H in 2HDMs, the exact decoupling limit of the extra scalar boson can no longer be 
taken, because of the unitarity constraints. Therefore, one can obtain an upper limit on the mass 
of the extra scalar boson [l^. The argument can be generalized to multi-Higgs doublet scenarios. 

The effects of a non-zero mixing can appear as deviations in the couplings of h with gauge bosons 
(hVV) and/or fermions {hff) from the SM predictions. Hence, the very fact that the latter can 
be accessed experimentally enables one to extract information on the underlying (pseudo)scalar 
dynamics above and beyond the mass scale directly probed at the collider concerned. Take the 
hVV couplings for example. They are presently (i.e., after LHC Run-I) constrained to be SM-like 
at the level of 10% or so. LHC Run-H (13 TeV) at standard integrated luminosity (300 fb“^) 


wou 


d improve this precision by about a factor of 2. Their expected accuracy will then be about 


2% [ill at the High-Luminosity (HL) LHC {l^ (with collision energy of 14 TeV and integrated 


luminosity of 3000 fb“^). Furthermore, the accuracy is expected to be about 0.4% {ll| at the 
International Linear Collider (ILC) with collision energy of 500 GeV and integrated luminosity of 
500 fb“^ (ILC500). Therefore, if a deviation in hVV couplings is discovered at future colliders, 
whenever this will occur, it will at the same time suggest the possibility of an enlarged Higgs sector 
and, if so, enable one to place an upper limit on the mass of extra (pseudo)scalar bosons. 

In this paper, we investigate constraints on dimensionless parameters in the scalar potential from 
the bounds imposed by perturbative unitarity and vacuum stability in 3-Higgs Doublet Models 
(3HDMs). The phenomenological importance of such extended Higgs scenarios, amongst many 
possible others, has been emphasized lately in the context of the latest LHC results and DM 
searches. For example, 3HDMs may shed light on the flavour problem, namely the unknown origin 
and nature of the three families of quarks and leptons, including neutrinos, and their pattern of 
masses, mixings and CP violation: it is possible that the three families of quarks and leptons 
could be described by the same symmetries that describe the three Higgs doublets {l3, 1^. Such 
a family symmetry could be spontaneously broken along with the EW symmetry, although some 
remnant subgroup could survive, thereby stabilizing a possible scalar DM candidate while for 
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certain symmetries it is possible to find a VEV alignment that respects the original symmetry of 
the potential which will then be responsible for the stabilization of the DM candidate |l^ . 

Although there are several types of 3HDMs depending on symmetries imposed in the scalar 


potential [16|]- having in view the aforementioned phenomenological requirements and a need 
for a simple and predictive formulation, recently, symmetric 3HDMs with one or two inert 
doublets have been investigated and found viable from both a DM and LHC perspective 


21 


2 ^. 


Furthermore, other than explaining current data, new regions of DM relic density and areas of 


parameter space yielding totally new LHC signals (with respect to, e.g., 2HDMs) open up [2ll-l23l|. 


Herein, we focus on the I(1+2)HDM and I(2+1)HDM, where the former [latter] model contains 
two[one] active doublets and one[two] inert doublets. We calculate the s-wave amplitude matrix for 
all possible 2-to-2 body scalar boson elastic scatterings in the high energy limit in both models, 
including contributions from the SM gauge bosons as the Nambu-Goldstone (NG) boson states 
due to the Equivalence Theorem Approximation (ETA)^. We then examine the allowed parameter 
space from the unitarity and vacuum stability (re-derived herein and found to be compliant with 


those in Refs. 


2ll |) bounds in addition to the constraints from the 5, T and U parameters. Then, 


as intimated, we proceed to apply the unitarity bounds to obtain an upper limit on the masses 
of A and H in the I(1-|-2)HDM in the case with a non-zero mixing between h and H, and study 
how the limit can be changed as compared to that in 2HDMs. We find that a stronger upper 
bound is obtained in the I(1-|-2)HDM than that in 2HDMs depending on the masses of the inert 
(pseudo)scalar bosons. Therefore, we can exclude the I(1-|-2)HDM when the masses of A or if will 
be found to be larger than the given upper limit for a fixed amount of the mixing between h and 
H. As for the I(2-|-1)HDM, since unitarity here is perfectly realized by the SM-like Higgs state h, 
we study what constraints are induced on the various parameters in the potential involving inert 
states) in order not to spoil the perturbative behavior of the only active scalar in this scenario. 

This paper is organized as follows. In Sec. H, we give the scalar potential and mass formulae 
for the scalar bosons. In Sec. HI, constraints from vacuum stability, perturbative unitarity and the 
S, T and U parameters are discussed. In Sec. IV, we delineate the allowed parameter regions from 
the constraints given in the preceding section. Finally, conclusions are given in Sec. V. 


^ In Ref. [2^, the similar approach has been applied to a 3HDM with a non-Abelian discrete S3 symmetry to obtain 
eigenvalues for the s wave amplitude matrix. Because the symmetry imposed in the Higgs potential in our paper is 
different from that in Ref. [20|], derived eigenvalues are also different. We have confirmed that our formula for the 
eigenvalues coincides with those presented in Ref. [2^ by taking appropriate replacements of the coupling constant 
in the potential. Detailed explanations are given in the end of Sec. III-B. 
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I(l-t2)HDM 

I(2-bl)HDM 

Doublet 


^2 

<Pi ‘^2 

^2 

- -b 

-b 

-b - - 

^2 

-b -b 

- 

-b -b - 


TABLE I: Charge assignment under the x Zi symmetry in the I(1+2)HDM and I(2+1)HDM. 


II. THREE HIGGS DOUBLET MODELS 

A. The scalar potential 


Among various classes of 3HDMs 
ments: 


In 


17l |. we discuss the one which satishes the following require- 


(i) it contains a SM-like Higgs boson; 

(ii) it contains a unique candidate of DM in a given mass spectrum; 

(hi) it does not have Flavour-Changing Neutral Current (FCNCs) at the tree level. 

In order to satisfy the requirements of (i) and (ii), we need at least one active and one inert 
doublet field. Thus, the (pseudo)scalar sector of our 3HDMs must be composed of one inert 
plus two active doublet fields; i.e., I(1-|-2)HDM or two inert plus one active doublet fields; i.e., 
I(2+1)HDM. To guarantee a doublet field being inert and the stability of a DM candidate, we 
require the corresponding VEV to be zero and impose an unbroken discrete symmetry onto the 
whole doublet sector. A DM candidate can then be obtained as the lightest neutral (pseudo)scalar 
component of the inert doublet field with the odd parity under Z 2 (the SM plus the Higgs/active 
(pseudo)scalar states being even). As a consequence, we also obtain that the inert (pseudo)scalar 
states can only appear in pairs in each vertex in the Lagrangian, so that they cannot couple to 
fermions directly. In the I(1-|-2)HDM, FCNC processes via neutral scalar boson mediation appear 
at the tree level similarly to the generic 2HDM case, because of the existence of two doublet fields 
with the same quantum numbers. The simplest way to avoid such FCNCs is imposing another 
discrete Z 2 symmetry Q which can be softly-broken in general (henceforth denoted as Z 2 ). Here, 
we consider both the I(1-|-2)HM and I(2-|-1)HDM with the two discrete symmetries Z 2 and Z-;^. 

® In the I( 2 -|- 1 )HDM, there is no FCNC at the tree level even when we do not impose the Z2 symmetry. However, 
by imposing the Zi symmetry throughout, the structure of the Higgs potential can be much simpler than that 
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To write down the scalar potential in the 3HDMs, let ^pi {i = 0, ...2) be isosipin doublet scalar 
fields, where the assignment of Z 2 x Z 2 charge is listed in Table [H The most general scalar potential 
is then given in the same form in both the I(1+2)HDM and I(2+1)HDM by 


n‘f(plpi + {nl2(p\p2+i^-C.) 

i=0,...2 
1 


+ 


Y + h.C.] 


i = 0,...2 

d, 


+ ai(iplip2)(PoPo) + CF2\P2‘fo\^ + + h.C.], 


+ + P2\PiP0?‘ + -j[p2,{P>\pof + h.C.] 


d, 


d, 


( 1 ) 


where p. 12 , A 5 , ps and <73 are complex parameters in general. Throughout the paper, we take 
these parameters to be real for simplicity, thereby avoiding explicit CP violation. We also ne¬ 
glect the possibility of complex VEVs, thus also removing spontaneous CP violation. Hence, the 
(pseudo)scalar helds can be parameterised as 


Pi = 


h: 


(f = 0,...2), 


( 2 ) 


where Vi are the VEVs of pi with the sum rule — (246 GeV)^. In the I(2-|-1)HDM, 

vi = V 2 = 0 and vq = v while in the I(1-|-2)HDM, vi, U 2 / 0 and vq = 0. 

The Yukawa Lagrangian in the I(2+1)HDM is given in the same form as that in the SM, where 
po gives all the SM fermion masses. In contrast, in the I(1-|-2)HDM, the structure of the Yukawa 
Lagrangian depends on the Z 2 charge assi gnm ents for the SM fermions. In general, there are four 


n 


independent choices of the assignment 3, 26] and these are the so-called Type-1 


ml27i. 


Vp^I, Type-X 


29, 3] 


strongly 


and Type-Y [271]. Various constraints from flavour 2^ and collider physics 
depend on this choice. However, the following discussion does not, though we will make sure to 
avoid phenomenologically dangerous regions of parameter space. 


B. Mass formulae in the I(1-|-2)HDM 


In the I(1-|-2)HDM, the inert sector is composed of one doublet. Therefore, the CP-even, CP- 
odd and pair of charged scalar states in the doublet field po correspond to the mass eigenstates. 


in the case without it, thereby rendering possible a semi-analytical treatment of unltarity violating processes. In 
this connection, it turns out that we can deal with the two models uniformly as seen in the scalar potential (given 
below in Eq. ©)• We thus consider the Z2 x Z2 symmetric potential in both 3 HDMs. 
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To avoid a confusion between the inert and active scalar states, we represent the inert states H^, 
Aq and Hq as and respectively. Their masses are calculated as 

= Mo + y [mi cos^ P + cTi sin^ /3] , (3) 

v‘2 

'^Vh ^ ^0 + y + P2 + Ps) cos^ P + {ai+a2 + as) sin^ /3] , (4) 

'^Va ^ y + M 2 - P3) cos^ /3 + (di + (72 - as) siy j3] . (5) 


The mass formulae for the active sector are completely the same as those in a generic 2HDM, 
so that we can directly apply the same mass formulae to the I(1+2)HDM. The mass eigenstates 
for the active scalar bosons are thus given as: 





with R{9) = 


( 6 ) 


cos 6 — sin 6 
sin 9 cos 9 

where and are the NG bosons which are absorbed by a longitudinal component of the 
and Z bosons, respectively, and h is the SM-like Higgs boson. The mixing angle /? is expressed by 
the ratio of the VEVs as tan/3 = ^ 2 /^ 1 . 

The squared masses of and A are then calculated as 


m^± = - y(A 4 + A 5 ), m\ = 


( 7 ) 


where 


= —4^- (8) 

sm p cos p 

The mass matrix for the neutral CP-even (scalar) states is expressed in terms of /3 and the usual 
mixing angle a as 


m\j = cos^(a — /3)M^^ -|- siy(Q; — /3)M|2 + sin2(a — /3)M^25 (9) 

m\ = siy(a — I3)MIi + cos^(a — /3)M|2 — sin2(Q; — /3)M^2i (10) 

9/rf2 

where the matrix elements are 

,;2 

MI^ = y (Ai cos"^ /3 + A 2 sin'^ /3) + y(A 3 + A 4 + A 5 ) sin^ 2/3, 

M 22 = sin^ /3 cos^ /3 [Ai + A 2 — 2 (A 3 + A 4 -|- A 5 )], 

2 2 

MI 2 = y sin2/3(—Ai cos^ /3 + A 2 sin^ /3) -b ^ sin 2/3 cos 2 / 3 (A 3 + A 4 + A 5 ). ( 12 ) 


7 




One of the most important variables in the active sector of the I(1+2)HDM, or equivalently in 
a 2HDM, is sin(/3 — a). This describes deviations in the SM-like Higgs boson h couplings to gauge 
bosons and fermions. Namely, the ratios of the hVV {V = Z) and hff couplings in the 

I(1+2)HDM to those in the SM are given by 


I(1+2)HDM 

9hVV 

9wv 


= sin(/3 — a), 


I(1+2)HDM 

3hff 


= sin(/3 — a) + Cf cos(/3 — a), 


(13) 


where the factors are cot /3 or — tan /3 depending on the type of fermion (i.e., up-type and down- 
type quarks and charged leptons) and the type of Yukawa interactions (i.e., the aforementioned 
2HDM Types)"^. The point here is that when we take the limit of sin(/3 — a) —i- 1, both the 
hVV and hff couplings become the same as those in the SM. We thus call this limit the SM-like 
configuration. Furthermore, by looking at Eq. m, we can see that only M 22 contains the term 
proportional to M^, while and M ^2 proportional to Thus, when we take the limit 
of m 2 —>■ 00 , the mixing angle given in Eq. dm) becomes tan 2 (a — /3) —>■ 0, which also gives 
m'jj —)• M 22 and —)• by choosing a — /? —>■ — 7 r/ 2 . Therefore, the extra Higgs bosons H^, A 

and H are decou pled due to the M^ term in their mass formulae, so that we can call this limit the 
decoupling limit [31 1. 


From the above discussion, we can reach an important conclusion regarding the nature of the 
Higgs sector in this model. The SM-like limit is naturally achieved when the decoupling limit is 
taken. In other words, the decoupling limit cannot be realized when we consider the case of a 
deviation from the SM-like limit, because we need a strong cancellation between and M 22 in 
the denominator of Eq. (ttH), in order to have tan 2{a — /3) 7 ^ 0. This would force the dimensionless 
coupling constants A* in the Higgs potential to acquire quite large values which must be excluded 
from the view point of ensuring perturbativity of the model. This in turn implies that there is an 
upper limit on the masses of the extra Higgs bosons when we consider the case of sin(/3 — a) 7 ^ 1. 
As mentioned previously, one of the aims of this paper is indeed deriving an upper limit for these 
in the I(1-|-2)HDM by using bounds from vacuum stability and perturbative unitarity, which we 
will discuss in the next section. 


^ The factor corresponds to 2T/with t/ being the third component of the isospin, where is given in Table II 

of Ref. [13]. 






C. Mass formulae in the I(2+1)HDM 


Next, we present the mass formulae in the I(2+1)HDM. In this model, (po corresponds to the 
SM Higgs doublet, so that and Aq correspond to the NG bosons and G^, respectively, 
while Hq is the SM-like Higgs boson h. 

The mass of h is given just like the SM Higgs boson as 

ml = 2\ov^. (14) 


In contrast, the inert sector is composed of two doublets, so that the masses for singly-charged, 
CP-odd and CP-even states are, respectively, given in the 2x2 matrix form. The mass matrix for 
the singly-charged scalar states is evaluated in the basis of H^) as 


Ml = 


pI + \pi 
fJ-12 


f^l2 


(15) 


P2 + 

and that for the CP-odd and CP-even states are given in the basis of {H^,H 2 ) and (^^^,^ 2 ), 
respectively. 


/ P2- P3 


|, -MH=-Mc + y| 

( P2 + P3 


V 0 

<72 — 0-3 ) 

1 2 \ 

{ 0 

0-2 + <73- / 


M\ = Mfj + 


The mass eigenstates can be defined by introducing three mixing angles, as follows: 


(16) 



= R{ec) 



= R{9a) 



= R{9h) 



(17) 


The squared mass eigenvalues for the singly-charged states ( 77 ^ , 772 ) expressed in terms of the 
elements of the mass matrix as 


= cos^ ec{Mc)u + sin^ 0ci-Mc)22 + sm29c{Mc)i2, 
m\ = sin^ 9c{Mc)u -h cos^ 6 'c(-Adc )22 - sin 20c(Me) 12 , 

H 2 


(18) 

(19) 


while those for the CP-odd and CP-even states ( 77 ^^, 77 ^^) and ), respectively, are obtained 

by replacing (0c,A4q) with {9a ,and {9h,M.‘\j)- Thus, the three mixing angles are given by 

In summary, the inert sector of the I(2-|-1)HDM can be described by nine parameters: i.e., six 
masses m ±, m„^ and 77 t,„ for i = 1,2 and three mixing angles 9x (for X = G,A,H). These 
nine physical parameters are expressed in terms of the nine parameters in the potential of Eq. m- 
i.e., pj, pI, PI 2 , Pi and Oi {i = 1,2,3). 


9 



III. CONSTRAINTS ON THE PARAMETER SPACE 


In this section, we discuss constraints on the parameter space in the I(1+2)HDM and the 
I(2+1)HDM. As theoretical constraints, we discuss bounds from vacuum stability and perturbative 
unitarity. As experimental constraints, other than ensuring sampling the 3HDM parameter spaces 
concerned away from experimental limits extracted via direct Higgs searches at LEP, Tevatron and 
LHC we also crucially take into account the so-called S, T and U parameters proposed by Peskin 
and Takeuchi 3^ in the high precision EW fits to LEP data 331]. 


For the discussion of vacuum stability and unitarity, there is no difference between in the 
I(1-|-2)HDM and I(2-|-1)HDM. In contrast, for the S, T and U parameters, different formulae are 
presented here in each of two models. 


A. Vacuum stability 

The Higgs potential should be bounded from below in any direction of the scalar boson space. 
The sufficient condition to guarantee such a positivity of the potential is given by the following 
requirements: 


Aq >0 , Ai > 0, A2 > 0, 

(21) 

\/ A1A2 -|- A3 -|- MIN(0, A4 -|- A5, A4 — A5) > 0, 

(22) 

V AqAi + Pi + MIN(0, p2 + ps, P2 — P 3 ) > 0, 

(23) 

V A0A2 -1- fJi -1- MIN(0, (72 + 0-3, (72 - ( 73 ) > 0. 

(24) 


B. Unitarity 


Constraints on the scalar self-coupling constants in the potential are extracted from the require¬ 
ment of S-matrix unitarity for all the elastic scatterings of two body scalar boson states. This idea 
to constrain a scalar self-coupling has been proposed in Ref. Q for the Higgs sector in the SM and 
it has been applied to constrain the mass of the SM Higgs boson. It can however be generalized to 
multi-doublet models, as well known 3^ . The requirement of the S'-matrix unitarity is translated 


into a relationship for the Jth partial wave amplitude aj. In the high energy limit, it can be 
expressed as 


Im(aj) = jajp. 


(25) 
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This implies that oj should be on a circle with radius and center of 1/2 and (0,1/2) in the complex 
plane, respectively. From the above relation, we can require the following condition for the tree 
level amplitude of aj as H 

|Re(aj)| < ( 26 ) 

In the high energy limit, each element of aj for all possible S 1 S 2 S' 3 <S '4 processes, where the 
Si's represent all (pseudo)scalar bosons in the model, is simply given by the (pseudo)scalar four 
point interaction. In that case, only the s-wave amplitude (J = 0) can contribute to the scattering 
process. We thus apply the inequality given in Eq. (|26|) to the case of J = 0. 

There are 30, 36 (= 18 x 2) and 12 (= 6 x 2) channels for electromagnetically neutral, singly- 
charged and doubly-charged states, respectively, where half of the singly- and doubly-charged 
states correspond to the charge conjugated states of the remaining half. Each channel classified 
by the difference in the Electro-Magnetic (EM) charge is orthogonal to the others, so that we can 
separately consider each of the channels. In order to distinguish these classes, we define o^, and 
as the s-wave amplitude for the neutral, singly-charged and doubly-charged states, respectively. 
In addition to the classification in terms of the EM charge, it is also helpful to distinguish the 
states by using the Z 2 x Z 2 quantum numbers of the intervening external states. Because we are 
now focusing on the (pseudo)scalar boson quartic interactions only, the softly-broken Z 2 parity can 
be considered as the exact quantum number. Under Z 2 x Z 2 , four types of (pseudo)scalar pairings 
are allowed: i.e., (^ 2 ,^ 2 ) = (-|-,-t-), (-|-,—), (—,+) and (—,—), where each state does not mix 
with any other one. We thus further decompose the matrices for the s-wave amplitude into four 
sub-matrices as follows: 



= diag (a++. 



(27) 

ao 

= diag (a++. 

a?-, 


(28) 

ao 

= diag (a++. 

a+-, 

a?_) . 

(29) 


Eor the singly- and doubly-charged states, we focus on the 18[6] channels with the positive charge, 
because the remained negative charged states are connected by the charge conjugation as explained. 

1. Neutral channels 

We first consider the neutral channels. There are 12 pairings, H^H^, AiAij^/2, HiHi/y/2 
and AiHi {i = 0, 1, 2) with (^ 2 ,^ 2 ) = (-|-,-|-). The 6 states with (-1 -,—) are H 2 , H 2 , 
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A 1 A 2 , H 1 H 2 , A 1 H 2 and H 1 A 2 . The other states with (—, +) and (—, —) are respectively given by 
replacing the subscript of the (+,—) states by X 1 Y 2 —)• XqYi and X 1 Y 2 —)• X 0 Y 2 . In fact, each 
sub-matrix in Qq can be made even simpler by choosing the following basis for the (-I-, -|-) states: 


^ , (30) 

^ (^H+H- - -^AA, - , (31) 

^ {AiAi - HiHi ), (32) 

T++)i+io = AiHi, with {i = 0, ...2). (33) 


For the (+, —) states, we can choose 

|T(^^_)i = i + H^H+ + ^ 1^2 + H 1 H 2 ) , (34) 

\^^.)2 = ^ + A 1 A 2 + Hi 112 ) , (35) 

|T(^_)3 = i (zH+H- - iH-H+ - A 1 H 2 + H 1 A 2 ) , (36) 

^ (-iH+H- + iH-H+ - A 1 H 2 + A 1 H 2 ) , (37) 

|T^^_)5 = ^(7li772 + 77i^2), (38) 

|T^^_)6 = ^(-^i^2 + 77i772). (39) 


Similarly, the basis for the (—,+) and (—,—) states denoted as |'I'(Y_,_)j and |'I'(Y_)j {j = 1...6) can 
be chosen by replacing the subscript in the above (-I-, —) states by X 1 Y 2 —)• XqYi and Xiy 2 —^ -^ 0 ^ 2 , 
respectively. In these basis, each sub-matrix in can be block-diagonal form with a 3 x 3 sub¬ 
matrix as the largest one. 


a 


N _ 
+-I- ~ 




a 


N 


X diag [(W|«) 1 , 


(40) 

Y^diag[n^_,n^_ 

4 


(41) 

Y^diag[ni+,n^+ 



(42) 

-^diag [n^ , 

4 

,n'^ ,n , 

,n^ , n® ] , 

(43) 


12 




where 


/ 


3Ao 2pi + p2 2 (Ti + (T 2 
2pi + p2 3Ai 2 A 3 + A 4 




y 2 ( 7 l+(T 2 2 A 3 + A 4 3 A 2 j 


I A, 


(iv: 


3x3a 

++ n 


0 P2 <72 
P 2 Ai A 4 
\a2 A 4 \2j 


( A, 


(iV3x3)3 ^ (^3X3)^ ^ 


0 P3 <73 
Ps Ai As 
\<73 As X 2 j 


_— A3 + 2A4 + 3 As, 

n\_ = A3 + As, 
n\_ = A3 — As, 

_= A3 + 2A4 — 3 As, 

n+_ = n%_ = A3 + A4, 
= pi+ 2p2 + 3p3, 
n?.+ = Pi +P3, 

n^-+ = Pi- P3, 
ni+ = Pi + 2p2 - 3p3, 
n^+ = n® + = pi + P2, 
ni_ = ai + 2(T2 + 3(73, 

n?._ = (7i +(73, 

= ( 7 i - (73, 

ni_ = ( 7 i + 2(72 - 3 ( 73 , 
n^_ = n® _ = (7i + (72- 


(44) 


(45) 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 


2. Singly-charged channels 

Second, we consider the 18 positive singly-charged channels. The sub-matrix has the 
form of 6 X 6 , while the sub-matrices a^_, a^_|_ and a(l_ have the form of 4 x 4. The basis 
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of can be expressed as {H^Ai, Hi) {i = 0 , ... 2 ), and that of a^_ can be expressed by 

A 2 ,H^H 2 , Ai, H 2 ). The basis for and a^_ are respectively given by replacing the 

subscript of the basis for o^_ by ( 1 , 2 ) —)■ ( 0 , 1 ) and ( 1 , 2 ) —>• ( 0 , 2 ). 

By choosing the following basis, we obtain a further simplified form of each sub-matrix; 


|4'5+). = T (,H+A + if+Hi) , 

(62) 

^ (ii7 +Ai - Hi) , with (z = 0, ...2), 

(63) 

i {iH+A2 + H+H 2 + iAiH+ + AiH+) , 

(64) 

|T^_)2 = i {iH+A2 + HtH2 - iAiH+ - A^H+) , 

(65) 

|T^_)3 = ^ {iH+A2 - HtH2 + iAiH+ - AiH+) , 

(66) 

|T^_)4 = ^ [iHtA2 - H+H 2 - iAiH+ + AiH+) . 

(67) 


The basis for the (—,+) and states denoted as and (k = can be 

chosen by replacing the subscript in the above (-I-, —) states by X 1 I 2 —^ -Toll and X 1 Y 2 —)■ -T 0 I 2 ) 
respectively. In this basis, each sub-matrix in Uq can be block-diagonal form as 

= Y^diag (C'^+^)2] , (68) 

a+_ = Y^diag [c \_, c\_ , ci^_, c|_] , (69) 

= Y^diag [c)_+, c?.+, ci+, ci+] , (70) 

a^_ = Y^diag [c)__, c?._, c^__ , ci_] , (71) 
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where 


(^3x^3; 

)i = {Nlf)2, {Clf)2 = {Nlfh, 

(72) 

c+- = 

nU 

(73) 

4- = 


(74) 

4- = 

n%-, 

(75) 

4- = 

A 3 — A 4 , 

(76) 

= 

n-+, 

(77) 

c4 = 

«-+> 

(78) 


n-+, 

(79) 

C-+ = 

Pi - P2, 

(80) 

ci_ = 

n--, 

(81) 

C2_ = 


(82) 

C^__ = 

n--, 

(83) 

ci_ = 

(7l — (72. 

(84) 


We note that c^_, and cl_ give the independent eigenvalues against those in the neutral 
states. 


3. Doubly-charged channels 


Finally, we consider the 6 positive doubly-charged channels. The sub-matrix has the form 
of 3 X 3 in the basis of {i = 0, ...2). Further, a^_, and a^_ have the 1x1 form 

for the H 2 , Hq and Hq H 2 states, respectively. In the doubly-charged channels, the above 
basis gives already the simplest structure of the matrix of the s-wave amplitude, written as 


o 

o 

a 

a 


D _ 
-i-i- ~ 

D _ 
— 

D _ 
-+ ~ 

D _ 
— 


1 

IbvT 

1 

Ibvr 

1 

Ibvr 

1 

Ibvr 


iDf+% 

d+-, 


h__ 


-d__. 


(85) 

( 86 ) 

(87) 

( 88 ) 
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They result in a 3 x 3 matrix with eigenvalues given as 


( 89 ) 

d_|_= n\_, d _ d _= n^_. (90) 


4 - Summary 


Consequently, we get 3 independent sub-matrices and 18 eigenvalues for qq. After we rename 
these, they are finally expressed by 



^ 3 Ao 2 pi-\-p 2 2 ai+o' 2 ^ 


^Ao P 2 < 72 ^ 


^Ao P 3 <73^ 

Xi = 

2 pi -|- p2 3 Ai 2A3 -|- A4 

, ^2 = 

P 2 Ai A4 

, ^3 = 

P 3 Ai As 


y2(7i -|- (72 2A3 -|- A4 3A2 j 


^(72 A4 A2y 


^(73 As X2) 


yf = ^3 + 2A4 ± 3A5, 

yf = Pi + 2p2 ± 3p3, 
y^ = ai+ 2fJ2 ± 3fJ3, 

2/J = A 3 ± As, 

yf = Pi± P 3 , 
y^ = ai± ( 73 , 

2/^ = A3 ± A4, 

y^ = Pi± p 2 , 

= (Tl ± (72- 


(91) 

(92) 

(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 
( 100 ) 


The condition in Eq. ()26l) is expressed in terms of the above variables as 


|xj| < 87r, (i = l,...9), (101) 

\yf\<87T, (i = l,...9), (102) 

where Xi are the eigenvalues of Xi, X 2 and X 3 . 

We note that all the eigenvalues of the matrices given in Eq. (fOT]l and those in Eqs. (lU^ - (llOOI) 
coincide with those given in Ref. , in which the s wave amplitude matrix has been derived in the 
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Ss symmetric version of the 3HDM, by taking the following replacement of the coupling constants 


Ai —)• 2 (Ai + A 3 ), A 2 —)• 2 (Ai + A 3 ), Aq —)• 2A8, 

A3 —)■ 2 (Ai — A3), A4 —)• 2 (—A2 + A3), As —)• 2(A2 + A3), 

Pi As, p2 —)• Ag, P3 2A7, 

< 7 i ^ As, (72 —)• Ag, (73 —)• 2A7, 

with taking A4 in Ref. Q to be zero. 


(103) 


C. S, T and U parameters 

We now discuss the S, T and U parameters in the context of our 3HDMs. The differences 
between the S, T and U parameters (S'shdm, T3HDM and Usubm) and those in SM (S'sm, Tsm and 
Usm), are defined by 


AS = 5'3hdm — 5'sm, 

AT = r3HDM — Tsm, 

AU = 173HDM — UsM- (104) 


Throughout this subsection, we use the shortened notations = sinX and = cosX. 

In the I(1+2)HDM, the differences are composed of two parts: i.e., the contribution from active 
scalar boson loops and that from inert scalar boson loops. The former one is the same as in 
the 2HDM while the latter one is the same as in the IDM, which is composed of one active and 
one inert doublet helds. Therefore, AS, AT and AU can be just obtained by summing the two 
contributions, active (^) and inert (/), as 


AS[I(1+2)HDM] = ASa + ASi, 
Ar[I(l+2)HDM] = ATa + ATi, 
At/[I(1+2)HDM] = AUa + AUi. 

The S, T and U parameters in the 2HDM and IDM have been calculated in Ref. 


(105) 


361] and in Ref. 


. fl. 
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respectively. Each of the differences is obtained as follows: 


ASa = Sp_^F’{m%]mH,mA) - F'{mz'^mH±,mH±) 




ATa = 


dvr 


+ c^_Q, F'{m%]mh, niA) + F'{m%; mu, mz) - F'{m%]mh, mz) 


+ 4m|c^_Q, G'{m%-,mH,mz) - G'(m|; m/i, mz) 
1 


lOvr^ 




F{{);mH±,rnA) + Sp_JyF{{);mH±,mH) - E(0;mA,m/f)] 


(106) 


+ c^_„ F{d]mH±,mh) + F{Q\mH,mw) + F{Q]mh,mz) 

- F(0; mh,mw) - E(0; mA,mh) - E(0; mH,mz) 

+ 4:G{0;mH,mw)+4:G{0;mh,mz) - 4:G{0;mh,mw) - 4:G{0;mH,mz) >, ( 107 ) 


AUa = {m^-,mH±,mA) - F {mz]mH±,mH±) 


+ s^-a[F'\miY-,mH±,rnH) - F'{mz'->mA,mH)] 


+ C/3_a 


F'{m^r;mH±,mh) + F'{m^;mw,mH) - F'{m^-,mw,rnh) 






— c 


fi—a 


F'{m%]'mA, ruh) + F'{mzl rnz,mH) - F'{mz]rnz, ruh) 


y /^2 




+ 4m^c^_Q G'{m‘^;mH, mw) - G'{m^-, mt, mw) 
G'{m%]mH.,mz) - G'{m%]mh, mz) 




ASi = 


ATi = 


AUj = 


47r - 


F {mz-,mn^,mr,J - F {mz;m^± ,m^±) 

V( 0 ;m^±,m^^) + F( 0 ;m^±,m^^) - F( 0 ;m^^,m^^) 


lOvr^aemU^ - 


1 

47r - 


F' (m^y; m^± , m ^^) + F' (m^y; m^±, ) 


- F'(m|;m^±,m^±) - F'{mz'^mn^,mr,^) 


( 108 ) 

( 109 ) 

( 110 ) 

( 111 ) 


where F'{my]mi,m2) = [F{my]mi,m2) — F(0; mi, m2)]/my and G'{my-,mi,m2) = 
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[G(my;mi,m 2 ) — G(0;mi,m 2 )]/my. The loop functions are given by 
F(p^;mi,m 2 ) = / dx (2x — l)(mi — m^) + (2x — 1) 

Jo L 


2 2 


in Ab, 


F(0; mi,m 2 ) = -(mi + m^) + 


2mJ:m% m 2 
in —, 


2' ^ ' m? — mo mi 


G{p 


'\rni,m2) = 

Jo 


= / dx in Ab , 


G(0;mi,m2) = ln(mim2) — 


mf + m^ m 2 


mj — m^ 


in-1, 

mi 


Ab = xmi + (1 — x)m 2 — x(l — x)p‘^. 


( 112 ) 

(113) 

(114) 

(115) 

(116) 


We note that the functions F and G given in the above are invariant by interchanging the second 
and the third argnments; i.e., F{jP‘-,mi,m 2 ) = F{p^;m 2 ,mi) and G(jP‘\mi,m 2 ) = G(p^;m 2 ,mi). 
We also note that ASi, ATi and AUi are obtained by taking sin(/3 — a) = 1 and making the 
replacement {mjj±,mA,mH) —>■ {m^±, m,y^,mr,^) in ASa, ATa and AUa- 

In the I(2+1)HDM, the differences purely come from the inert (pseudo)scalar boson loops. The 
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analytic expressions are given by 


^ 4^ I ~ ^ %± > %±) - P ^nt > %2±) 

+ ce^.eu ^'( 0 ; ) + ^'( 0 ; ) 


+ seA-e. 


AT = 


IGvr^ 




2 1 '^0C-0A 


+ Sec-dA 
+ ^9c-0h 


F'(0; ) - F'(0; ) 

F(0; m ^±, ) + F(0; m^±, ) 

F(0; m^±, ) - F(0; m^±, ) 

F(0; m^±, ) + F(0; m^±, ) 

+ [^( 0 ; %± > ) 

- ^dA-OH [-^(0; ^riA, ’ ) + ^(0; \ 

F(0; ) - F(0; ) 


— s. 


9a—Sh 


(117) 


(118) 


= —< - F {mz]rn^±,m^±) - F {mz]rn^±,rn^±) 


+ Ce^-dH [^ ) + F i'mw'^m^±, ) 

+ sec-9H [^' ; "ir,f > )-F'i^w-^ %± > ^Vh, ) 

+ C 0 ^_ 0 ^ [-F' (m^; m^±, ) + F' (m^; m ^±, ) 

+ S9^_eA %±, ^ ) 

- Cg^_g^ [F'{m%] rrin ^^, ) + F'(m|; ) 

9a—9h I F'(m|; ) - F'{m%\mn ^^, ) 


— s. 


(119) 


The deviations in t 


U = 0 are given by 


re S and T parameters from the SM predictions under rrih = 126 GeV and 


M 


AS = 0.05 ± 0.09, AT = 0.08 ± 0.07, 


( 120 ) 


where the correlation factor between AS and AT is +0.91. We use the above experimental values 
for the numerical analysis in the next section. 
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FIG. 1: Constraints on the I(1+2)HDM parameter space expressed in the the mn-mA plane induced by 
unitarity, vacuum stability and the S and T parameters in the case of mH± = mA, m^± = rririH, — 63 
GeV and sin(^ — a) = 0.98 (left panel) and 0.996 (right panel). We take to^± = 400, 500 and 550 GeV for 
the blue, green and red contours. We also show the result in the 2HDM as the black contour for comparison. 
The value of tan (3 is fixed to be 1 in the dotted contours while it is scanned over the range of 1 < tan /3 < 30 
in the solid contours. For all the plots, we scan the value of NP in the range of IVP = ± 1 TeV^. 

The outside regions from each contour are excluded by unitarity and vacuum stability. The light and dark 
shaded regions are also excluded by S and T parameters in the 2HDM and I(1+2)HDM, respectively. 



FIG. 2: Same as Fig.[T] but we take mfj± = mn instead of = ttia- 


IV. NUMERICAL ANALYSIS 

In this section, we discuss how the dimensionless quartic coupling constants in the scalar poten¬ 
tial can be constrained from perturbative unitarity, vacuum stability and the S and T parameters 
in our 3HDMs. Among these coupling constants, Aq and Ai ^,,,5 just describe four-point interactions 
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FIG. 3: Constraints on the I(1+2)HDM parameter space expressed in the ps-pi plane induced by unitarity 
and vacuum stability in the case of =mA=rnH{= m^± = 63 GeV and sin(/3 — a) = 

0.98 (left panel) and 0.996 (right panel). We take = 300[300], 500[500], 600[700] and 650[900] GeV for 
the black, blue, green and red contours in the left[right] panel. The value of tan/3 is fixed to be 1 in the 
dotted contours, while it is scanned over the range of 1 < tan^ < 30 in the solid contours. For all the 
plots, we scan the value of in the range of = rnjj ± 1 TeV^. Outside regions from the contours are 
excluded by unitarity and vacuum stability. 


among the inert (pseudo)scalar bosons in the I(1+2)HDM and I(2+1)HDM, respectively, and they 
do not affect to the active sector or the interaction between active and inert sectors at the tree 
level. We thus take these coupling constants to be zero in this section. We hrst show numerical 
results in the I(1+2)HDM and then in the I(2+1)HDM. 

A. I(1+2)HDM 

In the I(1+2)HDM, it is interesting to see how the constraints on the parameter space of the 
masses of active scalar bosons, A and H, can be modified by the requirement of unitarity and 
vacuum stability as compared to the 2HDM case. As we can see in the sections III-A and III-B, 
the additional coupling constants which do not enter in the 2HDM, such as pi and dj, appear in 
the conditions of vacuum stability and unitarity. Therefore, a non-zero value of these coupling 
constants can modify the constraints on the parameters describing the active (2HDM-like) sector 
i.e., A 1 -A 5 , which are translated into bounds on the masses of active (pseudo)scalar bosons. 

In the following, we fix the mass of either pA or pn to be = 63 GeV in order to satisfy 

the relic abundance of DM and take all the other masses of the inert (pseudo)scalar bosons 
to be larger than m/j/2. For definiteness, we choose p^ as the DM candidate, so that we take 
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fUri^ = m/i/2. To avoid a large contribution to the T parameter, we take m^± = tua or mij± = niH 
and m^± = For simplicity, we take cjj = pi {i = 1, ...3). In that case, the coupling constants 

Pi,...z are expressed in terms of the masses of inert scalar bosons as follows: 


^2 = 


1 


^3 = = -P2- 




( 121 ) 

( 122 ) 

(123) 


We can see that the sign of ps is now positive due to the assumption of = mh/2 (< nirf^)- 
The sign becomes negative when we consider the other case, of as the DM candidate, but the 
bounds from vacuum stability and unitarity do not depend on the sign of p^. From the vacuum 
stability condition of (1231) . the value of pi should be larger than 2|p2|) thus we take pi > 2|p2|- 
As we discussed in Sec. II-B, when sin(/3 — a) / 1 is taken, we cannot have the decoupling 
limit of the extra active Higgs bosons and thus we can obtain an upper limit of their masses. In 
addition, recall that the value of sin(/3 — a) describes the deviation in the hVV couplings as shown 
in Eq. ()13p at the tree level so that we may be able to find a relationship between the upper limit 
on the masses and the deviation in the hVV couplings. Because the l-a error of the measurement 




of hVV couplings is expected to be 2% at the HL-LHC and 0.4% at the ILC500 IllP . we take 
sin(/3 — a) = 0.98 and 0.996, respectively, as examples. 

In Fig. dl we show the constraint on the parameter space of the I(1+2)HDM mapped onto the 
fnH-fnA plane by unitarity, vacuum stability and the S and T parameters. We take sin(/3—a) = 0.98 
in the left panel and 0.996 in the right panel. The value of tan (3 is fixed to be 1 in the dotted 
contours, while it is scanned over the range 1 < tan/3 < 30 in the solid contours. For all the 
plots, we also scan the value of in the range of ± 1 TeV^. We check that these 

scanned ranges are enough wide to obtain the maximally allowed region on the parameter space 
from the unitarity and vacuum stability constraints. The blue, green and red contours show the 
cases of m^± = 400, 500 and 550 GeV, respectively. We also add a plot for the case in the 2HDM 
as the black contour for comparison. The regions outside each contour are excluded by unitarity 
and vacuum stability. In addition, the light and dark shaded regions are excluded by the S and T 
parameters in the 2HDM and I(1+2)HDM, respectively. We note that the contribution to AT is 
given to be zero due to mjj± = tua, so that the excluded region denoted by the light and dark area 
are from the constraint by AS". Besides, we conhrm that the difference in AS* due to the different 
choices of is negligible. 

In Fig. [21 we show the case of 'mjj± = mn with the rest of the setup being the same as in Fig. [T1 
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In this case, there is a non-zero contribution to AT which is proportional to cos^(/3 — a), so that 
the excluded regions denoted by light and dark shading are different from those in Fig. [TJ 

We find that the upper limits on niA and mn are related to each other and are getting stronger 
by taking larger values of m^±. We note that, if we take the limit m^±(= —)• the 

constraint from unitarity and vacuum stability approaches to that in the 2HDM. 

In Fig. [3l we show the constraints on the parameter space on the ps-pi plane by unitarity and 
vacuum stability in the cases of sin(/3 —a) = 0.98 in the left panel and sin(/3 —a) = 0.996 in the right 
panel, where the regions outside each contour are excluded. The values of tan /3 and are taken 
as done in Figs. [Hand[2j The black, blue, green and red contours show the cases of = 300[300, 
500[500], 600[700] and 650[900] GeV in the left[right] panel, where = mH± = uia = mH- We 
note that the parameter space shown in this plot is allowed from the S and T parameters. We 
can see that the region with pi < 2/93 is excluded by the vacuum stability condition in Eq. (|23p . 
Furthermore, the region above each contour is excluded by the unitarity constraints through the 
largest eigenvalue of the matrix Xi given in Eq. (j9ip . 

B. I(2-kl)HDM 

Next, we consider an application of our bounds from unitarity and vacuum stability in the 
I(2-|-1)HDM. As mentioned in the beginning of this section, we take A 1 -A 5 to be zero. In that case, 
we obtain the analytic expression of the eigenvalues for three sub-matrices of the s-wave amplitude 
given in Eq. (IMI) as 

of XiOi = diag ^0, ^ ± y^Ag -k (2pi -k ^ 2 )^ + (2fTi -k < 72 ) 2 ^ , (124) 

] J ^-^0 + P2 + ^2^ ) (125) 

O 3 X 3 O 3 = diag ^0, y ± y^Ag + pI + , (126) 

where Oi ,...3 are the orthogonal matrices to diagonalize Ai^,,, 3 . Now, Aq is determined by Aq = 
m|/(2u^) ~ -1-0.13 from Eq. (fTTp . so that the positive eigenvalues are important to get a stronger 
bound. 

In the following calculation, rjAi is taken to be the DM candidate with mass m/i/2. For simplicity, 
we consider a no-mixing case between two inert scalar states, which can be realized by setting 


O 2 X 2 O 2 = diag 
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FIG. 4: Constraints on the I(2+1)HDM parameter space expressed in the p^-pi plane induced by unitarity 
and vacuum stability. Outside regions from the contours are excluded by unitarity and vacuum stability. 

/^12 ~ 0^. In addition, to avoid the constraint from the T parameter, we take = rn^± and 

= rn^±. We note that results do not change if we replace both[either] and[or] 

with and[or] . In the above set up, we obtain the following relations for the quartic 

coupling constants: 


P2 = -\P3\, 0-2 = -|cr3|. (127) 

From the vacuum stability condition in Eqs. (j23p and (j24l) . we have 

Pi > 2|p 2|! <7i > 2|cr2|. (128) 

By the combination of the bounds from unitarity and vacuum stability, the case with pi > 2 \p2\ 
and cJi > 2 |cj2 | gives the largest allowed parameter space. In this setup, we find that the eigenvalue 
of Eq. (11241) gives the strongest bound. 

In Fig-U we show the allowed parameter space mapped onto the P 3 -pi plane in the I(2+1)HDM. 
Regions outside the contours are excluded by unitarity and vacuum stability. More precisely, the 
region below the contour with pi = 2 p^ is excluded by vacuum stability and the region above the 
contour is excluded by unitarity. 

® When we take pi 2 to be exactly zero, the lightest state among i Vh 2 is also stable. In order to have a 

decay mode of it, a non-zero value of /ri 2 is required. 
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V. CONCLUSIONS 


We have discussed the constraints on the parameter space of the I(1+2)HDM and I(2+1)HDM 
emerging from the bounds on vacuum stability, perturbative unitarity and the EW S, T and U 
parameters. To impose the unitarity bound, we have calculated the s-wave amplitude matrix 
for all possible 2-to-2 body (pseudo)scalar boson elastic scatterings in the high energy limit by 
means of the ETA, where the sub-matrices with 30 x 30, 36 x 36 and 12 x 12 forms for the 
electromagnetically neutral, singly-charged and doubly-charged channels, respectively, can be in 
block diagonal form with 3x3 sub-matrices as the largest component. We have first applied onr 
formulae for the unitarity and vacuum stability bounds to constrain the parameter space in the 
I(1-|-2)HDM, especially mapped onto the mn-ruA plane. We have found that larger excluded 
regions on snch a plane can be obtained in the case with sin(/3 — a) 7^ 1 as compared to those in 
2HDMs. We then have applied our procedure to the I(2-|-1)HDM on the ps-pi plane. Here, it has 
been clarified that the values ps > 8, 7, 6 and 2 are respectively excluded for any values of pi in 
the case of <72 = 0, —4, —6 and —8. 
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